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Подробно описывается алгоритм синтеза биортогональных условию кратномасштабного анализа. Приводятся примеры восстановления и 
сжатия сигналов на основе сконструированных вейвлетов. Описывается их сходство с базисом Карунена- Лоева. 


Введение 

Интенсивность и фаза оптического излучения, про- 
шедшего атмосферу со случайными неоднородностями, 
используются для извлечения информации об атмосфере 
или передачи информации через неё. В результате взаи- 
модействия с турбулентной средой распространения фаза 


и интенсивность волны становятся топологически слож- 
ными объектами. Происходит преобразование гладкого 
рельефа волнового фронта в прерывистую изломанную 
структуру со степенными особенностями и сингулярнос- 
тями. Для удобства и упрощения анализа эти сложные 
математические объекты, условно будем называть их 
случайными оптическими сигналами, обычно представ- 
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ляются в виде сумм ортогональных составляющих бес- 
численным количеством способов. И поскольку каждый 
раз система ортогональных функций, применяемая для 
разложения, бывает известна, то интенсивность или фаза 
оптической волны полностью определяются наборами 
весовых чисел для этих функций. Насколько адекватным 
реальности окажется восстановление сигнала, зависит от 
удачи выбора базиса. При решении задач восстановления 
поля волны возникает необходимость хранить большие 
массивы данных изображеніи, необходимость устранять 
шумы изображеніи, сократить время обработки данных 
и тем самым приблизиться к реальным временным мас- 
штабам физического процесса, отслеживать эволюцию 
частоты сигнала, вызванную различными временными и 
пространственными масштабами неоднородности. 

В данной работе мы не будем рассматривать поли- 
номиальные функции, гладкость которых заставляет их 
вести себя определённым образом, потому что сигналы 
с которыми приходится встречаться на практике хотя и 
непрерывные, но не дифференцируемые. Полиномиаль- 
ные функции не достаточно гибки, что бы отслеживать 
скачки и степенные особенности натурных сигналов. 
Поэтому в данной работе мы будем рассматривать 
вейвлет-базисы, обладающие фрактальным свойством 
самоподобия и позволяющие отслеживать сингулярности 
и изломы исследуемого сигнала. Замечательное свойс- 
тво самоподобия является причиной других полезных 
свойств, как то: свойство локальности, свойство быстрых 
преобразований, возможности высокой степени комп- 
рессии сигналов и изображений, возможность выделять 
сингулярные особенности и фрактальную структуру слу- 
чайного сигнала, возможность разложения исследуемого 
сигнала по масштабам неоднородности. 

Существует широкий класс вейвлет-базисов различ- 
ной природы, однако, их математическое описание слабо 
освещено в отечественной литературе. Количество отечес- 
твенной литературы по синтезу вейвлет преобразований не 
велико. Отметим некоторые из них. Полезной и интересной 
работой о возможностях приложения вейвлет преобразо- 
ваний является работа [1]. Наиболее информативными 
являются работы [2, 3] в которых приводятся классифи- 
кация вейвлетов и все необходимые теоремы по синтезу 
вейвлет преобразований. Недавно появился перевод книг по 
конструкции вейвлетов [4] и [5] . И наконец самым мощным 
источником информации о вейвлетах остаётся ІШетеР На 
сайтах ЬПр://\ѵ\ѵ\ѵ/тайі80й/сот/\ѵаѵе1еіЬйп1, ЬПр://р1ауіаіг. 
Яапіогс1.ес1іі/~\уаѵе1аЪ, \ѵ\ѵ\у.\ѵаѵе1е1.ог§, \ѵ\ѵ\ѵ.та{Ь.8рЬи. 
т/~с1тр/ можно получить информацию о программах, 
статьях, конференциях по теории и приложениям вейвле- 
тов. Последний из перечисленных является сайт Санкт- 
Петербургского семинара “Всплески и их приложение”. 
Для инженеров, интересующихся приложениями готовых 
программ по вейвлет-преобразованиям, можно порекомен- 
довать известные математические программные средства, 
как МАТЬАВ-6. 1, МаіЬетагіса-4 и МаіЬсасІ-200 1 (\Ѵаѵе1е1: 
ЕхШпзіоп). 

В данной работе мы проведём синтез новых вейв- 
лет-базисов и покажем на примерах их применение для 
обработки сигналов. Новые они потому, что ни один из 
приведённых ниже вейвлетов не относится к известным 


группам которые приводятся в существующей литера- 
туре по данной теме. Из широкого класса вейвлетов мы 
будем конструировать только биортогональные вейвлеты 
и только те, которые подчиняющиеся свойству кратно- 
масштабного анализа. Именно это свойство позволяет 
осуществлять эффективные быстрые алгоритмы разло- 
жения и восстановления. 

В отличие от синтеза ортогональных вейвлетов 
[6], синтез не ортогональных вейвлетов имеет гораздо 
большую свободу как по выбору формы, гладкости, 
симметрии, так и по таким важным критериям, как 
локализация и число нулевых моментов, то есть, класс 
не ортогональных вейвлетов гораздо шире. В этой 
работе мы подробно опишем алгоритм построения не 
ортогональных вейвлетов, которые имеют сходства с 
статистически оптимальным базисом Карунена-Лоева- 
Обухова [7]. 


Скейлинговые функции. Биортогональный вейвлет-базис 


При разложении сигнала Дх) в не ортогональный 
вейвлет-базис необходимо иметь два базиса, вейвлет 
Ф(х), по которому собствен"? раскладывается сигнал и 
двойственный ему вейвлет (х), с помощью которого 
определяются кс 11 ния, то есть 

ф (х)=Х>ЧѴ(Ч 

где Т-(Ъ=Г , '‘ч.(2-^-”) „ опкделяется 


как 


-со 

причём выполіу у ~т 

<ЧѴЛ',) Ѵ„8..„ _ (1) 

где 5 т . - символ Кронекера. Вейвлеты ^ и на- 
зываются биортогональными " образуются на основе 
скейлинговых функций ^ и ^ которые подчиняются 
двухмасштабным соотношениям (2) 

<рё)=ІУ<р(2*- к ), 

( 3 ) 


ф(Ъ=ІУф( 2 4 - к ). 

* =0 

№ 

Верхние индексы сѵмміггэвания N и являются 
носителями ф; м- ш ^ соответственно то есть 
К = 8иррф(Гп -зиррф. После определения коэф- 
(Ъициеунов и из формул (4) можно определить 
Ф и с помошью соотношений 

Ф(^)=І(-1У , У 1 Ф(2^- К ), (4) 

к =0 


^(Ъ=Х(-о +1, Ч,Ф( 2% - к ). 

к =о 

Таким образом, чтобы получить вешуг-ты ^ и 
необходимо определить коэффициенты и в раз- 
ложения (4). Для определения коэффициентов запишем 
разложения (А) к частотной области 

ф(к) = Ф(к/2)^у. ехр(- ^к/2), 
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$(к) = $(к / 2) і ]Г \ ехр(- ^к/2), 
2 к = 0 

Ф(к) = - 1 - } ФС ’ ) ехр(- * к) 5 * , 


2п 
где 

чения для фипътпов 


и введём обозна- 


1 


ЧД 


(к) = -2*- ехр(- <|Х к/2), (6) 

2 к = п 


^(к) = ^^ е хр(-^к/2). 

2 к = о 

лполнение условия ортогональности (3) для 'Р 
и * ’ эквивалентно выполнению следующего условия 
для филь™™ ги « 1 . 

Д С ) р (к) + д ( к + л) р (к + я) = 1, 

здесь черта означает знак комплексно^ """ляжения,. 

Если в качестве фильтра вьрГ~"" ' ” то решая 
уравнение (8) мы можем найти 

Прежде чем объяснить как решается уравнение (8), 
сделаем некоторое отступление, а затем вернёмся к его 
решению. 


Критерии синтеза вейвлет-базиса 




алгоритм синтеза вейвлета с требуемым числом нулевых 
моментов. 

Сконструировать вейвлеты 1 Р и ^ - значит найти 

й (к) ^(к). тз г ц (к) 

ѵ ’ и ѵ ' Выбираем ѵ ’ в виде 

ш(к) = С08”(к/2)Р( 8Іп 2 (к/2)), 

если п чётное, и в виде 

т(к) = ехр()к/2) со8 п (к/2)Р(8Іп 2 (к/2)), 

если п нечётное. " , іажении п - степень перво- 

С05 (— ) 

9 

го сомножителя ^ определяет число нулевых 
моментов в конструируемом вейвлете, а "'у'' некий 
полином. При выбрауь. ж "”"ыражении для будем 

искать решение для в виде 

( 12 ) 

(13) 


^ (к) = сов (к/2 ) 1 1 (зіп‘ і (к / 2)) 

если п чётное, и в виде 

и “(к) = схр( ф к/2)х 

ХС08 (к/ 2) П (8Ш : (к/2)), 

если п нечётное. 

Введём обозначеніи 

X 8ІП (к 2 ), 0-Ѵ 2,1 (к)=^(к), 
подставим (12) или (13) и соответственно (14) или (13) в 
(8), в результате получим уравнение Безу [8]: ( 14 ) 

В (Ъ = ”(' і ) П ?Х Ѵ = Г + ; )/2. п; 


п. 


( 8 ) 


В предыдущем разделе говорилось о том, что при 
сжатии изображения желательно иметь вейвлеты с 
большим числом нулевых моментов. Дадим краткое 
пояснение. * ;я в ряд 

фК 

если разложить сигнал 1'(х) в ряд Тейлора в окрестности 
точки к. то можно получить 

ф(^)=ф-( К )+ ф( К х^- К )+ 

1 

2 ѵ ! 

ЧД2 ^ -?сЬ 

Если мы умножим это выражение на ѵ 1^1 
и проинтегрируем то, очевидно, в разложении пропадут 
члены для I------- - чие 

\^(гт щ =о, = о,і, 2 ...м. 

В разложении (9) останутся высокочастотные 
члены при этом подавляются малоинформатив^іэд) 
низкочастотные тренды. При сжатии сигнала произ- 
водится сортировка (значимых) коэффициентов выше 
некоторого порогового значения. Очевидно, что, чем 
больше количество нулевых моментов тем лучше сжатие 
сигнала, то есть тем меньше останется коэффициентов 
для представления сигнала. В этом смысле симі^вд) 
ричный вейвлет-базис близок к оптимальному базису 
Карунена-Лоева [7]. Сходство проявляется в том, что 
базис Карунена-Лоева имеет некоррелированные, а вей- 
влет-базис имеет слабо коррелированные коэффициенты 
разложения, что является критерием для хорошего сжатия 
сигнала. Слабая коррелированность вейвлет-базиса объяс- 
няется тем, что носители базисных функций (замыкание 
множеств, на которых функции отличны от нуля) разне- 
сены в пространстве (или во времени). Возвращаясь к 
уравнению (6), скажем, что оно позволяет синтезировать 
вейвлеты с любым числом нулевых моментов. Приведём 


где 

урав 


Общее решение 


в (^) = 2 Х ( М + к -1,^Ч^Ф(^-1/2), 


(15) 


где С (К, к) - число сочетаний из п-элементов по к, а 
Дх) - любая нечётная функция. Первый и второй члены 
выражения (16), есть частное решение неоднородного 
уравнения и общее решение однородного уравнения 
Безу, соответственно. Выбирая различные N и Дх) 
при выбранном т( к ) в соответствии с (16^ "'щучим 
совершенно определённое выражение для * к *' При- 
ведём конкретный пример для чётного и нечётного N. 
Выбираем N=3 и Дх)=х+ х 3 , в этом случае имеем 

В(х)= 1+Зх+6х 2 +х 3 (х-1/2+(х-1/2) 3 ). (17) 

Определим корни выражения (18) и запишем его в 
факторна з 2 

Ѵ)= А П(* - 2К ^-)+ )’ 


где А - константа, которую можно опред( с ть, р с, 
делив с 9) на (18); в нашем случае А=1, =Ке( 

)+І I іп X ) ~ 3 пары комплексно сопряжённ'^” 

ней: '=-0,4763090739-і1,2303119751, 2 = ■» 
С С- 4 

3 = -0,2533865279 І0,3 024 8926 1^1 ' — 3 ’ 3 

= -1,47969571803 -і 1,224911 19812, «“ Теперь 
выбираем (18) 

Р(х)=( х 2 -2Ке(/ 2 )+|2 2 | 2 ), п=5, 
т( к )=ехр(]'к/2) соз 5 (к/2 )Р(8Іп 2 (к/2)), 


тогда (К ' П (Ъ = (^-2Ке(^)+|^ | 2 ) ХІ 
х(^-2Ке(^)+|^| г ), ° = 1, 


(19) 


^ (к) - ехр( ф к / 2) С08(к / 2 п (зіп 2 (к / 2)). 
Раскладываем ' ’ и ѵ ’ в ряды Фурье (7) и 
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Рис. 1. 
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получаем по десят: ’Д~Офиниентов и дляскс;- 
линговой функции и двойственной ей функции 
соответственно (рис. 1, а, б). Разложение в ряд Фурье 
можно осуществить, используя тригонометрические 
соотношения ѵ (20) 

( 21 ) 


СОЯ ( ? ) 


Я111 




І-* ^ 


( 22 ) 

4 нт ' (23) 

но гораздо проще разложение Фурье осуществлять в 

программно интегрируемой среде Майісасі. Далее, следуя 

алгоритму, описанному в [6], можно построить скейлинго- 

вые функции и соответствующие им вейвлеты. 


Графические зависимости для вейвлетов ^ и Ч' и 
их скейлинговых функций приведены на рис. 1, а, б. 
Для конструкции следующей пары вейвлетов следует 
в (19) принят т м " ~ 7, а^ П1 3 принять я = 3, при этом в 
разложениях * к ( и в ряд Фурье (7) мы пс *~ 
чим 12 коэффициентов и 10 коэффициентов 
(рис. 1, в, г). Соответствующие этим коэффициентам 
скейлинговые функции и вейвлеты приведены на рис. 
X в, г. Если принять п=5, я=1 и переместить корень 
1 из (21) в (19), то получим коэффициенты и соот- 
ветствующие им вейвлеты и скейлинговые функции, 
приведённые на рис. 1, д, е. 

Остановимся теперь на примере для чётного N = 4. 
В этом случае решение уравнения Безу будет иметь 
следующее решение 

В(х) = 1+4х+10х 2 +20х 3 +хД(х-1/2). 
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Рис. 2 


Выберем в качестве Р(х) нечётную функцию 
Р(х)= х+ х 3 . Находим корни В(х) и записываем в факто- 

риз° ІВ( , )=А(5 _^ ( , 2 _ 2Ке(Ч) + | Ч | 2) 

В этом случае получаем один вещественный корень 
х =-0,Гс "8926ГС '■ и 3 пары комплексно сопряжённых 
кор'^й = Ке( " )+]Іт( 2 .): ^ - 

1 = —0,9055499544 — і 1 ,58 1 306455, -> = з* 


; ’ = -8,1141156107-10- 2 -і0,372500б^ 

5= 1,9056374362 -іі, 5814099917, *“ 2 у ^ 
Теперь можно записать выражение (2ч) 

виде (12) и (14), соответственно, где * * и ( * 


есть 


различные сочетания сомножителей вьгр, ия,Т 

Раскладывая в ряды Ф' л “ьеё'д __ ътры к и ” 
получим коэффициенты “и " ,, , 

Возможные варианты решений приведеньг на 
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Рис. 3. Разложения сигнала по вейвлетам с разным числом коэффициентов разложения и с разным числом нулевых моментов вейвлетов 


рис. 2. Здесь мы привели только некоторые из возмож- 
ных решений для чётных и нечётных N. Отметим, что 
при фиксированном N класс возможных решений можно 
расширить, изменив вид функции Т(х), но при этом 
увеличиваются носители вейвлета и его скейлинговой 
функции. 

Рассмотрим теперь разложение сигналов по по- 
лученным базисам. На рис. 3 приведён модельный 
профиль фазы оптического сигнала (рис. 3, а) и его вей- 
влет-разложение (рис. 3, б-д). В качестве базиса были 
выбраны вейвлеты, приведённые на рис. 1, а - число 
нулевых моментов 8 = 3, и на рис. 2, д - число нулевых 
моментов, равное 8 = 6. На рис. 3, б, в приводится 
восстановление сигнала с компрессией для вейвлета 
с 8 = 6; рис. 3, б - 32 коэффициента из 256 и рис. 3, 
в - 128 коэффициента из 256. На рис. 3, г, д все то же 
для вейвлета с 8 = 3. Как и следовало ожидать, лучшее 


восстановление и лучшее сжатие сигнала (а, следова- 
тельно, и меньшее число коэффициентов в разложении) 
потребовалось для вейвлета с 8 = 6 с более высоким 
числом нулевых моментов. 
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